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Measurement of the chemical and isotopic composition of cosmic rays is essential for the precise understand- 
ing of their propagation in the galaxy. While the model parameters are mainly determined using the B/C ratio, 
the study of extended sets of ratios can provide stronger constraints on the propagation models. In this paper 
the relative abundances of the light nuclei lithium, beryllium, boron and carbon are presented. The secondary to 
primary ratios Li/C, Be/C and B/C have been measured in the kinetic energy range 0.35 — 45 GeV nucleon -1 . 
The isotopic ratio 7 Li/ 6 Li is also determined in the magnetic rigidity interval 2.5 — 6.3 GV. The secondary to 
secondary ratios Li/Be, Li/B and Be/B are also reported. These measurements are based on the data collected 
by the Alpha Magnetic Spectrometer AMS-01 during the STS-91 space shuttle flight in 1998 June. Our ex- 
perimental results are in substantial agreement with other measurements, where they exist. We describe our 
light-nuclei data with a diffusive-reacceleration model. A 10—15 % overproduction of Be is found in the model 
predictions and can be attributed to uncertainties in the production cross-section data. 
Subject headings: Cosmic rays - Acceleration of particles - 

Nuclear reactions, nucleosynthesis, abundances - Space vehicles 



1. INTRODUCTION 

The origin and properties of the charged cosmic rays (CRs) 
are one of the major subjects of modern astrophysics. Though 
experimental information comes from the analysis of the ar- 
riving fluxes, the understanding of the relation between ob- 
servational data and source properties requires a consistent 
picture of cosmic ray transport in the galaxy. Propagation 
calculations take into account the acceleration, energy losses, 
nuclear interactions, magnetic diffusion and convective trans- 
port of CRs through the galactic medium (Strong et al. 2007; 
Maurin et al. 2001). 

The propagation of CR nuclei is studied using the ratio of 
secondaries, which are created by fragmentation of heavier 
elements, to primaries, which are produced and accelerated in 
the astrophysical sources. The simple observation that the ob- 
served CR composition is different from that of rare solar sys- 
tem elements such as lithium, beryllium and boron proves the 
importance of propagation in the interstellar medium (ISM) 
in terms of fragmentation processes. In particular, the ratio 
B/C between boron and its main progenitor carbon is used to 
constrain quantities such as the average amount of interstellar 
matter traversed by CRs between creation and observation, 
or their characteristic escape time from our galaxy. In de- 
scriptions based on diffusion theory, the secondary to primary 
ratios are mainly sensitive to the energy dependence of the 
diffusion coefficient D. 

The light elements Li and Be are also of interest. Their 
abundances depend not only on interactions of the pri- 
mary species C, N and O, but also on tertiary contributions 
like Be— s>Li or B— >Li. Therefore the Li/C and Be/C ra- 
tios may provide further restrictions on propagation mod- 
els (De Nolfo et al. 2006). 

An accurate understanding of CR properties is also of im- 
portance in the search for exotic signals in rare components of 
the cosmic radiation, as the astrophysical background of any 
possible new physics signal must be estimated on the basis 
of the existing models (Salati et al. 2010; Moskalenko et al. 
2002). 
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The origin and evolution of the elements Li-Be-B is also 
a crossing point between different astrophysical fields: cos- 
mology, astroparticle physics and nuclear physics (Reeves 
1994). The model of Big Bang Nucleosynthesis (BBN) is 
able to produce only faint traces of nuclei up to A = 7. Al- 
most all the stars consume the relatively fragile Li-Be-B nu- 
clei in their thermonuclear core reactions. The galactic Li-Be- 
B are principally produced by the interaction of CRs with the 
ISM (Vangioni 2000). Despite many observations and pro- 
posed solutions there are still open questions about the high 
value of the 7 Li/ 6 Li ratio in meteorites and about the mea- 
sured overabundance of primordial lithium (Asplund et al. 
2006). 

The CR chemical composition has been extensively 
studied in a wide charge and energy range both on 
short duration balloon experiments (Webber et al. 1972; 
Orth et al. 1978; Lezniak & Webber 1978; Simon et al. 1980; 
Bukleyetal. 1994; Haryama & Shibata 2006), long du- 
ration balloon flights (Ahn et al. 2008; Panov et al. 2007) 
and space experiments (Webber et al. 2002; Engelmann et al. 
1990; S wordy et al. 1991). Isotopic composition measure- 
ments come mainly from space experiments such as the High 
Energy Telescopes on Voyager 1 and 2 (Webber et al. 2002) or 
the Cosmic Ray Isotope Spectrometer on the Advanced Com- 
position Explorer satellite (De Nolfo et al. 2003). Measure- 
ments in a wider energy range, but with poorer statistics, have 
been obtained by several balloon based magnetic spectrom- 
eters (Hams et al. 2004; Ahlen et al. 2000; Webber & Kish 
1979; Buffington et al. 1978). A new generation of balloon 
borne experiments is devoted to extending the measured range 
toward the knee, where the power law decrease in the flux 
of CRs appears to steepen, by means of transition radia- 
tion detectors and calorimetric techniques (Ave et al. 2009; 
Ahn et al. 2010; Panov et al. 2007). The current and upcom- 
ing space experiments PAMELA (Orsi et al. 2008) and the Al- 
pha Magnetic Spectrometer (AMS, Aguilar et al. 2002), rely 
on magnetic spectrometers, in which a magnetic field is used 
to bend the path of charged particles as they pass through. 
This allows accurate particle identification together with a 
precise momentum determination free from atmospheric in- 
duced backgrounds. 

In this paper we present the complete analysis of our mea- 
surement of the CR charge composition and energy spectra 
in the range of 0.35 — 45 GeV per nucleon. Our preliminary 
results were presented for the B/C ratio in Tomassetti et al. 
(2009a). The 7 Li to 6 Li isotopic ratio is measured in the 
rigidity region 2.5 GV to 6.3 GV (rigidity is defined as mo- 
mentum per unit charge, R = pc/Ze). These measurements 
take advantage of the large acceptance, accurate momentum 
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resolution and good particle identification capabilities of the 
AMS-01 spectrometer. 

2. THE AMS-01 EXPERIMENT 

The Alpha Magnetic Spectrometer (AMS) is a particle 
physics detector designed for the high precision and long 
duration measurement of cosmic rays in space. AMS-02 is 
scheduled to be installed on the International Space Station 
in 2010. The AMS-01 precursor experiment operated suc- 
cessfully during a 10 day flight on the space shuttle Discovery 
(STS-91). Data taking started on 1998 June 3. The orbital in- 
clination was 51.7° and the geodetic altitude ranged from 320 
to 390 km. 

The spectrometer was composed of a cylindrical permanent 
magnet, a silicon microstrip tracker, Time-Of-Flight (TOF) 
scintillator planes, an aerogel Cerenkov counter and antico- 
incidence counters (Aguilar et al. 2002). The magnet (inner 
diameter 1.1m) provided a central dipole field with an ana- 
lyzing power BL 2 = 0.14 Tm 2 . Six layers of double sided 
silicon microstrip tracker measured the trajectory of charged 
particles with an accuracy of 10 /im in the bending coordi- 
nate and 30 /im in the nonbending coordinate, as well as pro- 
viding multiple energy loss measurements. The TOF system 
had two orthogonal segmented planes at each end of the mag- 
net, and measured the particle transit time with an accuracy 
of ~90psec for Z > 1 ions. The TOF scintillators also pro- 
vided energy loss measurements up to \Z\ = 2. A layer of 
anti-coincidence scintillation counters lined the inner surface 
of the magnet. A thin carbon fiber layer was used as a shield 
to absorb low energy particles. 

A total of 100 million triggers were recorded during the 
mission. After the flight, the detector was calibrated with he- 
lium and carbon beams at GSI, Darmstadt, and with proton 
beams at the CERN PS, Geneva. The data of the precursor 
flight have provided significant results on galactic and trapped 
protons, electrons, positrons and helium nuclei (Aguilar et al. 
2002; Fiandrini et al. 2004, 2003). During the flight, nearly 
200,000 nuclei with charge Z > 2 were recorded by AMS- 
01. 

3. DATA ANALYSIS 

The physical characteristics of a detected CR particle are 
the arrival direction, the particle identity and its momentum. 
These quantities were reconstructed combining the indepen- 
dent measurements provided by the various detectors. The 
particle rigidity was obtained by the deflection of the particle 
trajectory. The velocity f3 = v/c was determined from the 
transit time between the TOF planes along the track length. 
The charge magnitude \Z\ was obtained by the analysis of 
multiple measurements of energy deposition (§3.2). 

The response of the detector was simulated using the AMS 
simulation program, based on GE ANT -3. 21 (Brun et al. 
1987) and interfaced with the hadronic package RQMD (Sorge 
1995). The effects of energy loss, multiple scattering, interac- 
tions and decays were included, as well as detector efficiency 
and resolution. Further details on reconstruction algorithms 
and detector performance are found in Aguilar et al. (2002) 
and references therein. 

During the STS-91 mission, the space shuttle Discovery 
docked with the MIR space station for about 4 days. The CR 
observations can be divided into four intervals: 

a) pre-docking, ~1 day of testing and checkout before 
docking with MIR; 
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FIG. 1. — Projection plot for downward going particles with measured 
charges Z > 2 when Discovery was docked with MIR. Note the relative 
deficit of events from a clearly defined region. Candidates originating from 
that region were removed by a geometric cut on the detector acceptance 
(dashed line). 



b) MIR-docked, ~4 days while the shuttle was docked 
with the MIR; 

c) post-docking, ^3.5 days pointing at fixed directions 

(V> = 0°, 20°, 45°); 

d) albedo, ~0.5 days pointing toward the Earth (ip = 
180°) . 

where ip is the angle between the AMS z-axis and the zenith. 
For the present work we considered data collected during in- 
tervals b) and c). We considered only particles traversing the 
detector top-down within a restricted acceptance of 32° from 
the positive z-axis. Data collected while AMS-01 was pass- 
ing in the region of the South Atlantic Anomaly (latitude: 5°- 
- 45° S, longitude: 5° - 85°W) were excluded. For data col- 
lected during b), a geometric cut on the detector acceptance 
was imposed, corresponding to the MIR's "shadow". In our 
previous work (Aguilar et al. 2005), this region was recog- 
nized through the secondary tt^ 1 and yU* produced by primary 
CRs interacting with the MIR modules. In the nuclear chan- 
nel (Z > 2), a flux deficit can be observed in the same region. 
This is shown in Fig. 1, where the region above AMS-01 is 
projected on to a X pr — Y pr plane using the standard trans- 
formations of the arrival direction X pr = — sin 6 cos <p and 
Y pr = sin 9 sin <p, where 8 and <j> are the polar and azimuthal 
angles in the AMS coordinate system of the incoming CR nu- 
clei. Approximately 85,000 Z > 2 candidates were selected 
with these criteria. 

3.1. Event Selection 

In order to reject events with large nuclear scattering or 
poorly reconstructed trajectories, a set of quality cuts was ap- 
plied to the selected candidates: 
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FIG. 2. — Top: normalized distributions of Q = x 2 1 f(R) ( a )> ^2 (b), and unassociated hit distance (c) for detected boron (solid circles) and carbon (open 
circles) nuclei. In the cases that no hits were detected in the silicon or of the track extrapolation falling outside of the active area, the distance is set in the 
d > 4 mm bin. The same distributions from the MC simulation data are superimposed (lines). Bottom: ratios between the corresponding distributions of B and 
C are shown for data (squares) and MC (gray lines). Cut thresholds are described by the arrows: events are accepted (a) to the left of the arrow, (b) between the 
two arrows, and (c) to the right of the arrow. 



1 . tracks fitted with large \ 2 were discarded according to 
a rigidity dependent requirement on Q = x 2 / f(R), 
where R is the measured rigidity; 

2. the particle rigidity was measured as Ri (i? 2 ) using 
the first (last) three hits along the reconstructed track. 
Tracks with poor agreement between R, Ri and R 2 
were removed according to cuts on A\ = (i?i — R) / R 
and A 2 = (R 2 - R) /R; 

3. consistency between j3 and R measurements was re- 
quired. This cut acted on the tails of the reconstructed 
mass distributions of the detected particles. 

4. silicon planes with no hit associated with the track were 
inspected in the region around the track extrapolation. 
The absence of additional hits was required within a 
distance d of 2 mm from the track extrapolation. 

As an indicator of mis-measured particles, a fraction of ~ 4 % 
of events with negative measured rigidity (R < 0) was found 
in the reconstructed data before applying the selection cuts. 
These events were used as a control sample to define the se- 
lection. The total selection efficiency was found to be sub- 
stantially charge independent over the whole energy range 
considered. In Fig. 2 the distributions of some quantities that 
were used for the selection are shown for boron and carbon in 
comparison with Monte Carlo simulation (MC). The agree- 
ment between data and MC in general is good; in the tails 
of the distributions, especially A 2 , the MC is less accurate in 
describing the measured behaviour (some few % of the sam- 
ple). However these features have similar magnitudes for the 
considered species and on average cancel out in the ratio, as 
shown in the lower panels of the figure. Similar conclusions 
can be drawn for the ratios Li/C and Be/C. 

3.2. Particle Identification 



Each elemental species is identified by its charge Z. The 
dynamical response of the TOF scintillators allowed particle 
discrimination up to Z = 2. The analysis of particle charge 
for Z > 2 events was performed by the study of the energy 
losses recorded in the silicon layers. The ionization energy 
generated by a charged particle in a silicon microstrip de- 
tector was collected as a cluster of adjacent strips. Tracker 
clusters were recognized online and then reprocessed with the 
reconstruction software. A multistep procedure of normaliza- 
tion of the cluster signals was performed (Tomassetti 2009b). 
The method accounted for saturation effects, electronics re- 
sponse, particle inclination and velocity dependence of the 
energy loss. The charge identification algorithm, applied to 
the corrected signals, was based on the maximum likelihood 
method which determines the most likely value of Z corre- 
sponding to the maximum value of the log-likelihood func- 
tion: 

fc=6 

L(Z) =log 10 {U Pl(x k ,l3)} (1) 
fc=i 

The fc-index refers to the six silicon layers, are the cor- 
rected cluster signals as observed in the fc-th layer. The nor- 
malized Pj (xk , /?) functions are the probabilities of a given 
charge Z with velocity (3 to produce a signal Xk in the fc- 
th layer. These probability density functions were estimated 
with a high purity reference sample of flight data, and de- 
scribe the energy loss distributions x k around a mean value 
Hz,/3 oc (Z 2 /f3 2 ) log(7) as expected by the Bethe-Bloch for- 
mula. This /3-Z dependence is shown in Fig. 3 for Z = 3 to 
Z = 8 (dashed lines) superimposed on the mean energy depo- 
sition from the measured tracker signals. Due to inefficiencies 
in charge collection, some energy losses may produce charge 
responses which do not carry reliable information on the par- 
ticle charge. For this reason, not all the six clusters were used 
to determine Z. For each reconstructed track, the most reli- 
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FIG. 3. — Mean energy loss in the silicon tracker as a function of the mea- 
sured velocity ft. The signal amplitude, shown in ADC count units, is the 
average signal of the best clusters {fci, &a, fes}. Different nuclear species fall 
into distinct charge bands. The fi-Z dependence of Hz, 8 i s superimposed for 
Z = 3 to Z = 8 (dashed lines). 
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FIG. 4. — Charge spectrum of the selected Z > 2 data. The signal ampli- 
tudes of Fig. 3 are equalized to (3 = 1 and shown in units of charge. Different 
nuclear species fall in distinct charge peaks. 



able set of three clusters = {^1,^2,^3} was selected by 
also using il as parameter of log likelihood funtion: 



L(z,n) = \og w {p h z ^p^-.p^} 



(2) 



where the free parameters to be determined are Z, ranging 
from 3 (lithium) to 8 (oxygen), and il, running over all the 3- 
fold combinations {ki, fc 2 , fc 3 } of the tracker signals. Though 
the final estimation of Z was provided by three tracker clus- 
ters, all the selected clusters were processed in Eq. 2. An iden- 
tification inefficiency of a few percent was achieved with this 
algorithm. On the basis of the final charge spectrum shown in 
Fig. 4, the determination of B and C was done with a charge 
resolution better than 0.14 charge units. 



3.3. Charge Contamination 

The charge assignment procedure was studied in terms of 
efficiency and contamination. Each nucleus of charge 3 < 
Z < 8 produces a charge estimation Z which can be related 
to its true impinging charge Z by using a set of coefficients 

The 6x6 matrix \\F\\ is diagonally-dominated, and each 
off-diagonal element F% represents the probability of a nu- 
clear species Z to be misidentified as Z due to interactions in 
the detector material and fluctuations of the energy loss: 



F z 



P(Z\Z) 



(3) 



Two different contributions produce a charge migration Z — > 
Z: 

A) After interacting in the upper TOF material, an incom- 
ing nucleus Z may fragment and physically turn into 
Z. These events were typically rejected by the antico- 
incidence veto. However, a fraction of them produces a 
clean track in the tracker, passing trigger and selection. 
Since Z < Z, the corresponding matrix is triangular. 

B) Fluctuations of the measured energy loss produce a 
nonzero Z —> Z migration probability. This charge 
spill over is typically symmetric and in most cases 
comes from adjacent charges, i.e. Z = Z ± 1. 

Both the effects were studied with Monte Carlo simulations. 
The second contribution was also estimated with the data. 
Possible charge migrations from/to the Z < 3 and Z > 8 
species were also considered, in particular to achieve a suit- 
able separation of the Z > 2 candidates from the back- 
ground of the more abundant helium flux. Helium-lithium 
separation was perfomed with the combination of both TOF 
and tracker energy depositions and studied with the help of 
data from the calibration beam test during which the detector 
was probed with pure helium beams at 0.75, 1.8, 3.6 and 8 
GeV nucleon -1 of kinetic energy (Aguilar et al. 2002). Con- 
tamination from He to the selected Z > 2 sample was esti- 
mated to be less than 10~ 4 of the helium sample. The effects 
of charge misassignment in the measured secondary to pri- 
mary ratios turned out to be smaller than the statistical uncer- 
tainties; this effect is included in the systematic errors (§4.3). 

4. FLUX DETERMINATION 

Before the secondary to primary ratios can be determined 
the spectra must be transformed into fluxes $. The differen- 
tial energy spectrum of the Z-charged particles measured by 
AMS-01 in the energy bin E of width AE is related to the 
measured counts AN Z by: 



§ Z {E) 



AN Z (E) 



A Z {E) ■ AT Z (E) ■ AE 



(4) 



where AT Z {E) is the effective exposure time (§4.1) and 
A Z (E) is the detector acceptance (§4.2). The relation be- 
tween the measured energy of detected particles and their true 
energy were studied using unfolding techniques (D'Agostini 
1995). Making use of an appropriate energy binning, the fi- 
nal smearing effects turned out to be barely appreciable and 
very similar for the variuous species, i.e., the overall effect is 
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FIG. 5. — Rigidity resolution of the AMS-01 silicon tracker estimated with 
measured data from the flight and the calibration test beam with carbon nuclei 
in comparison with the MC simulation. 



mostly canceled out in the bin-to-bin ratios. The rigidity res- 
olution curve of AMS-01 is shown in Fig. 5 for carbon nuclei 
measured during the flight and in the beam test, together with 
the simulation results. 

4. 1 . Exposure Time 

By the term exposure time AT Z (E), we mean the effective 
data taking duration with AMS-01 capable of receiving a new 
trigger from a galactic particle of charge Z in the energy in- 
terval AE. The detection is influenced by the DAQ livetime 
fraction a and the geomagnetic field modulation. Though no 
atmospheric or trapped particles are expected in our Z > 2 
fluxes, the position dependent geomagnetic cutoff introduces 
different distortions of the measured energy spectrum for dif- 
ferent nuclear species. Hence only particles with rigidity 
greater than R Th = 1.25 • R c were accepted, where Rc is 
the rigidity cutoff obtained from the Stormer estimation in the 
corrected geomagnetic coordinates for 1998 (Smart & Shea 
2005). The entire orbit was divided into 45,000 time intervals 
of 8t ~ 10 s. For each time bin 5tk, the mean values of R^' 1 
and «fc (trigger livetime fraction) were computed. The expo- 
sure time AT Z (E) was then calculated for each particle Z, 
rigidity R and energy E e AE as: 



AT- 



Oik 



(R,Rl h ) 



(5) 



where R = (|) y/E 2 + 2M P E, H = 1 for R > R Th , and 
H = for R < R Th , i.e., the sum is restricted to only the 
time intervals above cutoff. The result describes the effective 
exposure of AMS-01 to cosmic rays in that energy bin com- 
ing from outside the magnetosphere. Note that AT Z {E) as 
a function of the kinetic energy per nucleon is different for 
different A/Z isotopes; since A/Z is slightly different for the 
elements under consideration, their corresponding exposure 
times do not completely cancel in the flux ratios. 

4.2. Acceptance Estimation 



The detector acceptance A z (E) of Eq. 4 is the convolution 
of geometrical factors with the energy dependent efficiency 
e, including trigger, reconstruction and selection efficiencies. 
The acceptances of each species A z were determined with the 
Monte Carlo simulation (Sullivan 1971). 2-10 8 trajectories 
were generated in the energy range 0.2 — 50 GeV nucleon -1 , 
according to an isotropic distribution and a power law mo- 
mentum distribution oc p _1 . The recorded MC events were 
sent through the same analysis chain as for the measured data. 
Background fluxes were also simulated, up to Z = 8. 

4.3. Systematic Errors 

Though systematic uncertainties arising from many steps 
of the analysis are suppressed in the ratios, differences in the 
trigger efficiencies of the two species are present, as expected, 
from the charge dependence of delta ray production and frag- 
mentation effects in the detector material. The spill over from 
adjacent charges also produces net effects on the measure- 
ments. 

The MC determination of the acceptance gave a statistical 
uncertainty of ^1—4%, increasing with energy. A Z (E) de- 
creases with increasing energy and charge, due to the trig- 
ger conditions against interacting particles. In the events col- 
lected from the flight, the dedicated unbiased dataset 9 did not 
give enough statistics for an energy dependent validation of 
the Z > 2 trigger efficiency with measured data, hence this 
estimation must rely on simulation results. Since CR inter- 
actions with the detector material can play an important role 
in the observed behaviour of the trigger response, a GEANT3- 
based MC simulation does not guarantee a model independent 
description of the experimental setup. Therefore, this study 
was performed with three different models. A Vi rtual MC 
application (Hrivnacova et al. 2003) was developed for AMS- 
01. The simulation natively supports the particle transport 
codes GEANT3 and GEANT4 (Agostinelli et al. 2003). An 
additional interface with FLUKA (Fasso et al. 2005) was de- 
veloped (Oliva 2007; Tomassetti 2009b). Since the quanti- 
ties that enter in to our measurements are the ratios rjx be- 
tween the corresponding acceptances A Ll /A c , A Be /A c and 
A B /A c , the systematic errors on the trigger efficiency were 
estimated observing the scatter between the three different 
models. Results of this approach are summarized in Fig. 6. 
From the figure the charge dependent behaviour of the trigger 
response is clear; the energy dependence of the trigger effi- 
ciency is more pronounced for higher charges. Good agree- 
ment was found between GEANT4 and FLUKA; nevertheless, 
for each ratio, the full envelope of the three model curves was 
considered as the uncertainty band of the trigger efficiency. 

Another source of uncertainty comes from the isotopic 
composition of the measured nuclei. In order to compute 
bin-to-bin ratios, all the elemental fluxes were determined 
in a common grid of kinetic energy per nucleon E. The 
rigidity to energy conversion required the A/Z ratio of each 
charged species. To first approximation A/Z » 2 for 
the light Z > 2 nuclei under study; however, according 
to existing measurements and model calculations, more re- 
alistic isotopic mixtures can be considered for the Li-Be- 
B group (Hams et al. 2004; Ahlen et al. 2000; Webber et al. 
2002; Buffington et al. 1978). Hence we assumed isotopic 
fractions 6 Li/( 6 Li+ 7 Li) = 0.5 ± 0.15, 7 Be/( 7 Be+ 9 Be) = 
0.65±0.15, and 10 B/( 10 B+ n B) = 0.35±0.15, accounting for 

9 One out of 1000 events was recorded with the AMS-01 unbiased trigger. 
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these uncertainties in the estimation of the systematic error. In 
the left panels of Fig. 7 we cross checked our measured ratios 
with independent measurements obtained with the TOF sys- 
tem, which provided a direct measurement of the velocity and 
hence kinetic energy per nucleon. Results from the tracker 
measurement and the TOF measurement are consistent. The 
right panels of Fig. 7 show the reconstructed mass distribu- 
tions for the Li-Be-B events above cutoff with /3 < 0.95. 
Though the AMS-01 mass resolution did not allow an event- 
to-event isotopic separation, the data are in good agreement 
with the MC simulation which contains the assumed mixtures. 

A detailed summary of the uncertainties on the measured 
ratios Li/C, Be/C and B/C is presented in terms of relative 
errors in Fig. 8. The sources of uncertainties taken into ac- 
count are statistics of collected events, MC statistics, contam- 
ination due to charge mis-identification (§3.3), uncertainty on 
isotopic composition and trigger efficiency uncertainty. 

5. RESULTS AND DISCUSSION 

5.1. The Li/ 3 Li Isotopic Ratio 

A composition fit technique for determining the isotopic ra- 
tio 7 Li/ 6 Li in four rigidity intervals between 2.5 and 6.3 GV 
was performed (Zhou 2009). The best fits were obtained 
by leaving the amplitudes for both isotopes free in all four 
rigidity intervals as shown in Fig. 9. The fits gave unique 
minima for the 7 Li/ 6 Li ratios in each rigidity region, pro- 
viding an average ratio of 1.07 ± 0.16 in the region 2.5- 
—6.3 GV. The results are shown in Fig. 10 with the previous 
experimental data (De Nolfo et al. 2003; Webber et al. 2002; 
Hams et al. 2004; Ahlen et al. 2000; Webber & Kish 1979; 
Buffingtonet al. 1978); the latter are converted into units of 
rigidity, allowing a direct comparison with our data. Point 
values are listed in Table 1. Our measurement agrees with the 
previous data and extends them to higher energies. In sum- 
mary, our new measurements of 7 Li/ 6 Li extends to 6.3 GV of 
rigidity (^2 GeV nucleon -1 of kinetic energy) with a con- 
stant ratio of about equal abundance of both isotopes. 

5.2. Secondary to Primary Ratios 

Results for the Li/C, Be/C and B/C ratios are pre- 
sented in Fig. 11 with the existing experimental data 
(Webber et al. 1972; Orthetal. 1978; Lezniak & Webber 
1978; Engelmann et al. 1990; Ahn et al. 2008; De Nolfo et al. 



TABLE 1 
The isotopic ratio 7 Li/ 6 Li. 



Rigidity (GV) 2.5-3.2 3.2-4.0 4.0-5.0 5.0-6.3 

7 Li/ 6 Li 1.11 ±0.16 1.02 ±0.13 1.10 ±0.16 1.03 ±0.18 
X 2 /d.f. 35.1/35 46.1/39 41.2/43 42.6/43 



2003). The energy range 0.35 — 45 GeV nucleon -1 is limited 
by selection inefficiencies below ~0.35GeV nucleon -1 and 
lack of statistics above 45 GeV nucleon -1 . The error bars in 
the figure represent the sum in quadrature of statistical errors 
with the systematic uncertainties. These results are summa- 
rized in Table 2. 

Our B/C ratio measurement agrees well with the results 
from the first flight of CREAM in 2004 (Ahn et al. 2008) 
and with the data collected by HEAO-3-C2 (Engelmann et al. 
1990) from October 1979 and June 1980. The Be/C ratio is 
consistent, within errors, with the HEAO data, but not with 
balloon data Orth et al. (1978). Our Li/C data have unprece- 
dented accuracy in a poorly explored energy region. In com- 
parison with balloon data from Orth et al. (1978), our data in- 
dicate a quite different trend in the high energy part of the 
Li/C ratio. 

In these ratios, the main progenitors of boron nuclei are pri- 
mary cosmic rays (CNO). On the contrary, the abundances 
of Li and Be depend also on secondary progenitors Be and B 
through tertiary contributions like B— >Be, Be— !>Li and B— >Li. 

However, the observed shapes of the measured ratios Li/C 
and Be/C suggest their suitability in constraining the propa- 
gation parameters: their decreasing behaviour with increasing 
energy is a direct consequence of the magnetic diffusion ex- 
perienced by their progenitors, while the characteristic peak 
around ~1 GeV nucleon -1 is a strong indicator of low energy 
phenomena like stochastic reacceleration or convective trans- 
port with the galactic wind. 

5.3. Secondary to Secondary Ratios 

For sake of completeness, in Fig. 12 and Table 4 we re- 
port the secondary to secondary ratios Li/Be, Li/B and Be/B. 
Though their values can be derived from the secondary to pri- 
mary ratios of Table 2, a dedicated analysis of uncertainty was 
done for these channels. These quantities are less sensitive to 
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TABLE 2 

The secondary to primary ratios of Li/C, Be/C and B/C between kinetic energies of 0.35 

AND 45 GEV/NUCLEON. THE COLUMNS INDICATE: ENERGY INTERVALS E 1 -E 2 , REFERENCE VALUES 
{E)", MEASURED RATIOS, STATISTICAL ERRORS AND SYSTEMATIC ERRORS. 



AE 


<E> 


Li/C 


STAT 


SYS 


Be/C 


STAT 


SYS 


B/C 


STAT 


SYS 


0.35 


-0.57 


0.45 


0.187 


±0.009 


±0.007 


0.094 


±0.007 


±0.005 


0.323 


±0.013 


±0.005 


0.57 


-0.92 


0.73 


0.200 


±0.007 


±0.008 


0.101 


±0.005 


±0.005 


0.334 


±0.010 


±0.006 


0.92 


-1.50 


1.18 


0.203 


±0.006 


±0.008 


0.106 


±0.004 


±0.004 


0.339 


±0.008 


±0.008 


1.50 


-2.44 


1.92 


0.184 


±0.004 


±0.009 


0.105 


±0.004 


±0.006 


0.302 


±0.007 


±0.011 


2.44 


- 3.97 


3.11 


0.159 


±0.004 


±0.010 


0.103 


±0.004 


±0.007 


0.265 


±0.008 


±0.012 


3.97 


-6.45 


5.06 


0.139 


±0.005 


±0.009 


0.091 


±0.005 


±0.007 


0.234 


±0.010 


±0.012 


6.45 


- 10.5 


8.22 


0.108 


±0.006 


±0.009 


0.082 


±0.006 


±0.007 


0.199 


±0.012 


±0.015 


10.5 


- 17.0 


13.36 


0.090 


±0.008 


±0.007 


0.073 


±0.009 


±0.006 


0.175 


±0.019 


±0.013 


17.0 


-27.7 


21.72 


0.083 


±0.014 


±0.007 


0.078 


±0.015 


±0.008 


0.141 


±0.030 


±0.011 


27.7 


-45.0 


35.30 


0.088 


±0.026 


±0.007 


0.064 


±0.023 


±0.010 


0.126 


±0.043 


±0.013 



a The reference energy (E) is computed as the geometric mean yJE\ E2 . 



the propagation parameters, because all the secondaries have 
the same origin and undergo similar astrophysical processes. 
These ratios are also found to be less influenced by solar mod- 
ulation. Hence, these relative abundances maximize the ef- 
fects of the nuclear aspects of the CR propagation: decays, 
fragmentations and catastrophic losses. As pointed out in 
Webber & Soutol (1998), the Be/B ratio is sensitive to the ra- 
dioactive decay of the beryllium (numerator) which enriches 



the boron flux (denominator); thus a precise measurement of 
the Be/B ratio provides constraints on the galactic halo size as 
well as the 10 Be/ 9 Be ratio. The observed increasing behaviour 
of the Be/B ratio is due to the relativistic lifetime dilation of 
the unstable Be. However, the migration 10 Be^ 10 B involves 
only a small fraction of the elemental flux, and the present 
data are still affected by sizeable errors. 
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Fig. 8.— 



Relative errors on the measured ratios Li/C, Be/C and B/C. 
The contributions from data statistics, Monte Carlo statistics, contamination 
on the charge identification, isotopic composition and trigger efficiency are 
shown as a function of the kinetic energy per nucleoli. For the sake of graph- 
ical clarity, points are connected with lines. 



5.4. Propagation Calculations 

To describe our results, we made use of GALPROP 
v5 0.1p 10 , a numerical code that incorporates all 
the astrophysical inputs of the CR galactic trans- 
port (Strong & Moskalenko 1998). We described the 
propagation in the framework of the diffusive-reacceleration 
model, that has been very successful in the description of the 
cosmic ray nuclei fluxes (Moskalenko et al. 2002). GALPROP 
solves the diffusion transport equation for a given source dis- 
tribution and boundary conditions for all the galactic CRs. In 
this model, the magnetic diffusion is described in terms of a 

rigidity dependent diffusion coefficient D = (3D (R/R ) S . 
The magnitude of the diffusion, related to the level of 
hydromagnetic turbulence, is controlled by the parameter 
Do which fixes the normalization of D(R) at the reference 
rigidity i?o- The spectral index 5 is linked to the density of 
the magnetic irregularities at small scales. Reacceleration 
is an energy gain of charged particles due to scattering on 
hydromagnetic waves moving at Alfven speed V a in the ISM. 
This process is described as a diffusion in the momentum 

10 http : / /galprop . stanf ord. edu 



TABLE 3 
GALPROP SETTINGS. 



Parameter 


Name 


Value 


diffusion - normalization 


D 


5.85-10 28 cm 2 s- 1 


diffusion - index 


5 


1/3 


reference rigidity 


Ro 


4GV 


Alfven velocity 


V a 


30kms _1 


injection - index 


a 


2.43 


halo height 


11 


4 kpc 


solar modulation 


4> 


450 MV 



space and controlled by the parameter V a . The source spectra 
q(R) are assumed to be have a pure power law dependence 
in rigidity, i.e. q(R) oc R~ a . The boundary conditions 
are expressed by imposing free escape at the galactic halo 
boundaries, in particular by the halo height H. Energy 
losses and catastrophic losses over the interstellar medium 
are included. The code makes use of an up to date nuclear 
reaction network, including decay rates and fragmentation 
cross sections for all the relevant channels. Semiempirical 
models for cross section calculations are tuned with mea- 
sured data, where available (Moskalenko & Mashnik 2003). 
Equilibrium solutions are provided for the local interstellar 
spectra (LIS) of all the galactic CRs. In our description, we 
used a 2D cylindrically symmetric model of the galaxy, For 
the heliospheric propagated fluxes, we adopted the force field 
description of the solar wind (Gleeson & Axford 1968). We 
used <fi = 450 MV as modulation parameter, consistent with 
1998 June (Wiedenbeck et al. 2005). The parameter values 
of our GALPROP settings are listed in Table 3. 

In Fig. 11 and 12 our data are compared with the model 
calculation in the heliosphere (solid lines). The LIS ratios 
(dashed lines) are shown for reference. Our results for the 
B/C ratio, the Li/C ratio and the 7 Li/ 6 Li isotopic fraction of 
Fig. 10 are described quite well within the uncertainties. It 
is difficult, however, to accomodate the Li and B description 
with the Be/C ratio by only means of astrophysical parame- 
ters. Beryllium appears to be overproduced in the model by 
a factor ^10—15 %. This discrepancy is also apparent in the 
previous measurements. The ratios Be/B and Li/Be of Fig. 12 
also indicate this feature, whereas the Li/B ratio is well de- 
scribed by our GALPROP tune. 

It is worth noting that the production of Li and Be involves 
very complex reaction chains, due to the multichannel charac- 
ter of the CR fragmentation. Moreover, cross sections for the 
tertiary processes LiBeB + ISM — > LiBe are less well known 
than for the interactions CNO + ISM ->• LiBeB. Spallation 
contributions like 11 B— > 9 Be are measured at only a few en- 
ergy points. Other channels like 10 B— >- 9 Be, 14,15 N— > 9 Be or 
1011 B— > 7 Be are not measured at all and rely on extrapolated 
parametrizations (Moskalenko & Mashnik 2003). Spallation 
processes involving interstellar helium are even less well un- 
derstood. The lack of cross section measurements limits the 
model predictions to uncertainties of ^10-20% in the Li/C 
and Be/C ratios (De Nolfo et al. 2006). Understanding frag- 
mentation is a key factor in establishing final conclusions con- 
cerning cosmic ray propagation. 

6. CONCLUSIONS 

We have presented a new measurement of light CR 
nuclei composition in the energy range from 0.35 — 
45GeV nucleon -1 with the AMS-01 experiment. The iso- 
topic ratio 7 Li/ 6 Li has been measured between 2.5 and 6.3 
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FIG. 9. — Lithium mass distributions and best fit composition in four rigidity regions. The amplitudes of the simulated distribution of 7 Li and 6 Li are free 
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TABLE 4 

Secondary to secondary ratios of Li/Be, Li/B and Be/B between kinetic energies of 0.35 
and 45 gev/nucleon. the columns indicate: energy intervals e1-e2, reference values 
{e)", measured ratios, statistical errors and systematic errors. 



AE 


<E> 


Li/Be 


STAT 


SYS 


Li/B 


STAT 


SYS 


Be/B 


STAT 


SYS 


0.35- 


0.57 


0.45 


1.935 


±0.157 


±0.225 


0.577 


±0.030 


±0.015 


0.291 


±0.023 


±0.014 


0.57- 


0.92 


0.73 


1.979 


±0.1 18 


±0.227 


0.598 


±0.024 


±0.014 


0.300 


±0.017 


±0.012 


0.92- 


1.50 


1.18 


1.889 


±0.085 


±0.221 


0.599 


±0.020 


±0.017 


0.311 


±0.013 


±0.012 


1.50- 


2.44 


1.92 


1.747 


±0.066 


±0.216 


0.610 


±0.020 


±0.027 


0.348 


±0.013 


±0.019 


2.44- 


3.97 


3.11 


1.544 


±0.065 


±0.202 


0.603 


±0.024 


±0.035 


0.384 


±0.015 


±0.026 


3.97- 


6.45 


5.06 


1.528 


±0.084 


±0.203 


0.594 


±0.032 


±0.037 


0.387 


±0.022 


±0.027 


6.45- 


10.5 


8.22 


1.310 


±0.105 


±0.183 


0.539 


±0.045 


±0.042 


0.404 


±0.033 


±0.033 


10.5- 


17.0 


13.36 


1.223 


±0.157 


±0.155 


0.515 


±0.069 


±0.036 


0.416 


±0.057 


±0.022 


17.0- 


27.7 


21.72 


1.062 


±0.238 


±0.135 


0.588 


±0.151 


±0.045 


0.553 


±0.139 


±0.023 


27.7- 


45.0 


35.30 


1.352 


±0.442 


±0.218 


0.699 


±0.280 


±0.056 


0.514 


±0.202 


±0.060 



il The reference energy (E) is computed as the geometric mean \/E\Ei- 



GV of rigidity. This work is aimed at investigating the origin 
and the physical properties of the galactic cosmic rays. 

The study of high charged ions with AMS-01 required the 
development of an improved charge identification algorithm. 
A thorough analysis has been made to understand the CR in- 
teractions in the spectrometer, its instrumental response and 
the orbital environment. Our results for the B/C ratio agree 
well with data collected from the HEAO-3-C2 experiment and 
with the more recent measurement from CREAM. The Li/C 
ratio has been measured with unprecedent accuracy. A 10- 
— 15 % overproduction of Be is found in the model predic- 
tions, that describe well Li and B. This is consistent with the 
lack of cross section measurements, that limits the model pre- 
dictions up to uncertainties of ^20 % for light-nuclei. 

The astrophysical interest of light nuclei cosmic ray fluxes 
measurements, both for propagation studies and in the search 



for exotic phenomenon, has been well established in the CR 
community. We expect, with AMS-02, to measure them with 
high precision over wide energy ranges in the near future. 
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FIG. 1 1 . — Results for the secondary to primary ratios Li/C, Be/C and B/C from this work (solid circles) and previous measurements (Webber et al. 1972; 
Orthetal. 1978; Lezniak & Webber 1978; Engelmann et al. 1990; Ahn et al. 2008; De Nolfo et al. 2003). Model calculations are also reported, from the 
GALPROP diffusive-reacceleration model predictions for interstellar (LIS) and solar modulated (</>=450 MV) cosmic ray fluxes (Strong & Moskalenko 1998). 
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FIG. 12. — Secondary to secondary ratios Be/B, Li/Be and Li/B from this work (solid circles) and previous measurements (Orth et al. 1978; Lezniak & Webber 
1978; Engelmann et al. 1990). Model calculations are also reported, from the GALPROP diffusive-reacceleration model predictions for interstellar (LIS) and solar 
modulated (<^=450 MV) cosmic ray fluxes (Strong & Moskalenko 1998). 



